Entanglement localization by a single defect in a spin chain 
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We discuss the effect of a single diagonal defect on both the static and dynamical properties of 
entanglement in a spin chain. We show that entanglement localizes at the defect and discuss its 
localization length, arguing that this can be used as a mean to store entanglement. We also show 
that the impurity site can behave as an entanglement mirror and characterize the bouncing process 
in terms of reflection and transmission coefficients. 
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Entanglement generation and distribution is a prob- 
lem of central importance in performing quantum- 
information tasks, like teleportation and quantum cryp- 
tography. As entanglement between parties is cre- 
ated by means of a direct interaction, methods are re- 
quired to transfer either the entangled particles or their 
state. Indeed, the problem of designing quantum net- 
works enabling efficient high-fidelity transfer of quan- 
tum states has been addressed by a number of authors, 
IH 00 [100 especially focusing on the requirement of 
minimal control. This means that state transfer should 
be possible without performing many control operations 
such as switching on and off the interactions, measuring, 
or applying encoding and decoding procedures. 

In this respect, spin systems provide ideal models to 
study entanglement properties as they are naturally em- 
ployed as qubit registers and exploited as quantum chan- 
nels (or coherent data bus). Spin chains with fixed in- 
teractions have been considered and solid state im- 
plementations have been already proposed [2j. It has 
also been shown that perfect transfer can be achieved 
by performing local measurements by using several 
spin chains in parallel 0, by employing local memories 
at the receiver side 0, or by means of a spatial modu- 
lation of the spin coupling strengths |fj. For this latter 
case, the effect of imperfections (static errors) in the en- 
ineering of the spin system has been analyzed in Ref. 
where a scaling of the transmission fidelity has been 
obtained in terms of the degree of imperfection. This is 
an example in which the disorder has been taken into ac- 
count. In fact, considerable attention has been devoted 
to disordered spin systems. A threshold in the coupling 
has been found for the onset of entanglement between a 
bulk impurity and its neighborhood [8(. Also it has been 
shown that a faster (super-ballistic) distribution of entan- 
glement between a central node and a series of distributed 
sites occurs if the sender resides in a disordered region . 
Furthermore, perturbativc, numerical, and Bcthe-ansatz- 
based investigations of entanglement between two defects 
has been performed for a disordered anti-ferromagnetic 
spin- 1/2 chain with anisotropic exchange coupling, (lfij . 
A possibility of tuning the ground state entanglement 



by a single off-diagonal impurity in the anisotropic XY 
model has also been considered jTll |. 

Diagonal disorder has also been considered in the quan- 
tum information context, but only the wave function lo- 
calization issue has been addressed in terms of the inverse 
participation ratio 12], while entanglement localization 
has been completely overlooked, apart from some consid- 
erations reported in Ref. 0] . 

It is well known that disorder can lead to a spatial 
localization of the electronic wave function and that, 
in particular, this occurs whenever impurities exist in a 
tight binding model. In fact, a single impurity suffices 
and, in this paper, we discuss how this translates 
into an entanglement localization. 

We consider the simplest example of a system display- 
ing localization; namely, a spin chain with a single di- 
agonal defect, that is, a chain subject to a local static 
magnetic field which is equal at every spin location ex- 
cept for a single site, where a field defect is present. We 
analyze both the static and dynamical properties of such 
a system, and, in particular, we discuss the case in which 
the ground state of the chain is a localized state, showing 
a very peculiar spatial distribution of quantum correla- 
tions. We also discuss how this localized state affects the 
possibility of sending entanglement through such a 'de- 
fected' channel. In particular, we show that a mirror-like 
effect occurs at the defect during entanglement propa- 
gation. This achieves a (partial) control of the entangle- 
ment dynamics, as, starting from a homogeneous system, 
the modification of the single-qubit level spacing permits 
the generation of an "artificial defect" . 

The model we study is a 1-D isotropic XY spin-i closed 
chain, placed in an external magnetic field which is ho- 
mogeneous everywhere but for a single (defect) site I 
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where h is the magnetic field, acting on every site but for 
site I, where the field is h + e. Also, J is a ferromagnetic 
coupling constant, while a a are the Pauli matrices. We 
will consider the impurity term as a perturbation and 
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take as a reference the results obtained for the homoge- 
neous XY Hamiltonian [II 03. Since [i?,E*<4] = °> 
the total Hilbert space can be divided into invariant sub- 
spaces labelled by the number of spins which point up- 
wards. We will concentrate on the and 1-excitation 
subspaces and denote by |n) the state with the upward 
spin at site n. 

Thorough, we consider the case 2h > J. This implies 
that the unperturbed problem (e = 0) has the factorized 
ground state lO)®^ 4 " 1 , whose energy we re-scale to zero. 
In the Hilbert space sector with one reversed spin, the 
unperturbed energy basis is 
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In this sector, the Hamiltonian JJJ is equivalent to a 
tight binding problem with one diagonal impurity. This 
can be solved with the help of the Green's function tech- 
nique (lif. 

The unperturbed Green operator is 



G (z) = 
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In the thermodynamic limit {N — * oo), the index k be- 
comes continuous and the states \k) give rise to a contin- 
uous energy band, with Eu <G [2h — J, 2h + J] = lb- In 
this continuum limit, the matrix elements of Go between 
localized states can be written 



G (r,s;z) = ± ' , z ^ I b (4) 
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where a; 

The Green operator for the full hamiltonian can be 
expressed in terms of the unperturbed one, and, due to 
the simple form of the interaction hamiltonian, the D yso n 
series can be re-summed giving an exact expression 
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G(z) = G (z) + G (z)\l) 
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The knowledge of the Green function ||BJ| allows us to ob- 
tain the spectrum of the hamiltonian It turns out 
that the energy band of the continuous spectrum is unaf- 
fected, as for E E [2h — J,2h + J] the Green function has 
a cut in the complex plane. The associated eigenstates 
are given by \V{E)} = J2n a n( E )\ n ), with 



ln(E) = 



1 



a i\e\\n-l\ 



i \ smf 



, (7) 



where a = ^ and cos6 = (2h — E)/J. 
Besides a distortion of the states within the band, 
the perturbation produces the appearance of a local- 
ized eigenstate, whose energy is given by the pole of 
the Green function JJjJ, determined by the equation 
2cGo(l,l;Ei oc ) = 1. This state lies above or below the 
band, depending on wether a is greater than zero or not, 
and its energy is Ei oc — 2h =p Jyl + a 2 . Explicitly, the 
state is given by |*z oc ) = J2 n & « l n )' witn 
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where £(a) = — ln (Vl + a 2 — \a\). is the spatial 

extension (or localization length) of |J[jpc) around the 
impurity site I (£ is obtained in Ref. |l8| for the XXZ 
model). For j-J==g < 1 with a < (that is, for negative 
and large enough a), Ei oc is lower than zero and |^j oc ) 
becomes the ground state of the system. 

We analyze the entanglement content of this localized 
state. In particular, we concentrate on the pairwise quan- 
tum correlations shared by the qubits at sites i and j as 
measured by the concurrence 

La- it is easy to show 
that the concurrence between these sites is given by the 
residue of the associated matrix element of the Green op- 
erator (to simplify the notation, we set / = and measure 
distances along the chain from the impurity site) 
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As a result of the presence of the magnetic field defect, 
pairwise entanglement is only present for spins residing 
within a distance l/£ from the impurity and it decays 
exponentially with the separation. In this sense, we state 
that entanglement is localized in the neighborhood of the 
defect. As the localization length becomes less than the 
site spacing for a < (1 — e 2 )/2e ~ —1.175, entanglement 
is highly localized even for small values of the defect field. 

We now analyze the transmission of a quantum state 
along the chain. In particular, we are interested in char- 
acterizing entanglement distribution; namely, the possi- 
bility of using the chain to send one partner of a max- 
imally entangled pair. A related problem is the trans- 
mission of a single qubit state along the chain. To send 
one qubit, we imagine that at site s, the sender prepares 
a generic single qubit state |^ s ) which is aimed to be 
retrieved at site r (where, in general, the mixed state 
p r (t) is extracted at time t). A measure of the quality of 
the chain as a transmission channel is given by fidelity, 
F(t) = (tp s \ p r (t) \ip s ), averaged over all possible messages 
\ip s ). On the other hand, to send entanglement, the spin 
s is initially prepared in a singlet state with an external 
(un-coupled) qubit. The entanglement transmission can 
be characterized by the (final) concurrence between the 
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FIG. 1: (up-left) Localization length as a function of \a\. 
(down-left and right) Ground state concurrence between the 
impurity and the spin at site j. Notice the exponential de- 
crease as a function of the distance. 



external qubit and the one residing at the receiving site 
r, denoted by C r (t). 

In fact, the two quantities F and C are related 0, 
as both can be expressed in terms of the transmission 
amplitude to send an excitation from s to r, f rs (t) = 



-iHt 



s) . Explicitly, the averaged fidelity is (F) (t) 



(\frs(t) + 1| 2 + 2)/6, whereas the concurrence is given by 
C r (t) — \frs\- In our case, the transition amplitude is 
given by a complicated combination of matrix elements 
of the Green functions, 
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where, within the integral, E = 2h— J cos 9, while 



(E) 



G^\i,j;E) 



,±i\B\\i-j\ 



1 -qG^' (1,1; E) ±ism\e\-a' 



The amplitude f rs given in eq. I|1U|) contains two terms. 
One contribution comes from the localized state; it has 
the same structure of the concurrence C rs given in eq. 
(0. The other contribution is written in terms of an in- 
tegral over the branch cut region of the Green function, 
describing the evolution of the states within the contin- 
uous band. While the former goes to zero with the am- 
plitude of the defect a, the latter does not, as it reduces 
to the un-pcrturbed transition amplitude when a — ► 0. 
This is given by the Bessel function of order r — s and 
argument r = Jt,J r - s (r), see Ref. [l^. The main effect 
of a finite defect field in the second term is a distortion of 



the amplitude, with the appearance of contributions de- 
scribing scattering at the defect during propagation from 
s to r. 

In the limit of a > 1 , a series of simple remarks can be 
drawn. First of all, the localization length goes to zero, 
and therefore the localized state is really concentrated on 
the impurity site. As a result, if the qubit residing there 
is initially entangled with an external one, the entangle- 
ment doesn't propagate at all, and the concurrence has 
the simple form (for s = I = Q): 
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Secondly, the transmission through the defect is very 
poor. Indeed, when the sender is on one side of the im- 
purity (located at I = 0), while the receiving site is on 
the other one, at a distance d — |r| + |s|, the concurrence 
is given by 



Ct 
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(12) 



which decays both with the distance between the sites 
and with time as r~ 3 / 2 . 

Finally, the defect can reflect entanglement. Indeed, 
when both the sender and the receiver are on the same 
side (say s,r > with r > s), the concurrence is given 
by the combination of three terms: 
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This relation has a simple interpretation: the entangle- 
ment can either propagate directly from s to r (a process 
weighted by the un-perturbed transition amplitude given 
by the Bessel function J r _ s ), or it can propagate to the 
defect and then be reflected to the site r (thus travelling 
a distance 2s + r — s = r + s). This latter process has an 
amplitude obtained by i) taking the Bessel function J r + S 
which describes propagation along this distance, and 2) 
subtracting (with its own relative phase) the amplitude 
for entanglement transmission through the defect. 

Although not easily seen from the full expression of 
eq. (|10|l . these two main features of entanglement local- 
ization and reflection at the impurity site, are present 
also for moderate values of a. This can be seen in Fig. 
(0), where the case a = — 2 is illustrated. By increasing 
a, these effects become more and more pronounced and 
the secondary propagation lines which are present in Fig. 
(j2J tend to disappear. 

To better characterize the effect of the impurity on 
entanglement propagation, we can introduce reflection 
and transmission coefficients. Indeed, we can exploit 
the fact that states with one reversed spin saturate the 
CKW inequality [2(|, which in our case implies that 
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FIG. 2: Time evolution of the concurrence C r for a = —2 and 
two different initial conditions. The left plot shows the case in 
which the sender resides at the defect site (I — s = 0), whereas 
the right plot shows the case s = —5. In the first case, the 
entanglement with the external qubit stays localized at the 
sender site, while in the second one, entanglement propagates 
almost symmetrically from the sender site up to the defect, 
where it is reflected back. 



dynamics of the chain. Indeed, the qubit at the defect 
site can either maintain entanglement for long time or 
work as an entanglement mirror. These two properties, 
obtained by modifying the local magnetic field of just one 
qubit, can be exploited to achieve a control on the prop- 
agation of pairwise quantum correlations in distributed 
qubit systems. Indeed, we have shown that entanglement 
can be stored for long times at the defect site and that 
its propagation can be modified so that it can be directed 
towards the receiver. 



This work has been done within the framework of the 
PRIN2005029421 project. 
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FIG. 3: Transmission and reflection coefficient as defined by 
Eq. (I14L displayed as a function of a. 



J2 r Cr (t) = lj and use the tangle as a probability dis- 
tribution to define (to be specific, we consider the case of 
s < 0) 

T= lim £c*(i), R= ]imJ2c 2 r (t). (14) 

r>0 r<0 

The transmission and reflection coefficients are shown in 
Fig. as functions of a. 

Entanglement reflection at the defect is completely ac- 
counted for by the presence of the localized state. Indeed, 
since the overlap between the initial state and l^joc) is ex- 
ponentially small, entanglement propagates almost freely 
up to the defect site, but then the spin-reversal excitation 
cannot enter there. As the Hamiltonian describes nearest 
neighbor interaction, subsequent spins cannot receive the 
excitation and are, therefore, almost excluded from the 
dynamics. As a result, entanglement is reflected back- 
wards. 

In summary, we discussed the effect of a single mag- 
netic field defect on entanglement in a spin chain. Wc 
found that this impurity produces a state with local- 
ized entanglement, whose presence strongly affects the 
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